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The  adsorption  of  the  chemical  warfare  agent  simulant  dimethyl  methylphosphonate  (DMMP)  and  the  real 
agents  Sarin  and  VX  on  the  y-AhOi  surface  has  been  studied  using  density  functional  theory.  The  focus  is 
primarily  on  two  different  environmental  effects,  namely,  surface  hydroxylation  and  photoexcitation  due  to 
terrestrial  solar  radiation.  Cluster  models  for  the  hydroxylated  surface  have  been  examined  in  detail,  focusing 
on  the  chemical  and  electronic  structure.  The  energy  for  formation  via  dissociative  adsorption  of  H20,  the 
density  of  states  of  the  occupied  cluster  orbitals  and  the  OH  deprotonation  energies  have  been  compared  with 
results  from  two-dimensionally  periodic  slab  calculations  and,  where  available,  with  experimental  data.  For 
all  three  species,  adsorption  on  an  OH-free  surface  occurs  via  an  Al(Td)~ 0=P  dative  bond  to  an  unsaturated 
tetrahedral  A1  (Td)  site.  For  the  hydroxylated  surface,  OH  sites  which  are  3-fold  coordinated  to  A1  are  more 
reactive  than  one-fold  coordinated  sites,  in  agreement  with  experiment.  In  hydrogen-bond  formation,  the 
phosphonyl  O  atom  is  favored  over  other  active  centers  (e.g.,  an  alkoxy  O  atom);  however,  dative  bonding 
remains  the  most  stable  mode  of  adsorption  when  OH-free  Al(7j/)  sites  are  available.  The  tertiary  amine 
group  in  VX  is  found  to  be  sufficiently  basic  to  extract  an  H+  from  an  acidic  OH  site.  Even  on  a  fully-OH- 
covered  surface,  only  a  single  OH— 0=P  bond  is  formed,  together  with  weaker  CH— O  bonds  involving 
alkyl  groups.  These  have  been  analyzed  using  the  Atoms  in  Molecules  theory.  The  threshold  for  electronic 
excitation  of  DMMP  or  Sarin,  either  in  the  gas  phase  or  adsorbed  on  y-Al2O  i,  is  found  to  occur  in  the  vacuum- 
ultraviolet,  well  beyond  the  upper  limit  of  the  terrestrial  solar  spectrum  (~4.5  eV).  For  VX,  on  the  other 
hand,  the  threshold  is  at  ~4.2  eV  in  the  gas  phase  and  shifts  slightly  to  the  red  when  adsorbed.  Adsorption- 
induced  shifts  in  the  threshold  transition  energies  have  been  analyzed  in  terms  of  those  in  the  initial  and  final 
states. 


1.  Introduction 

The  application  of  quantum-chemical  methods  to  the  study 
of  the  phosphonyl-based  chemical  warfare  agents  (CWAs)  has 
attracted  a  significant  amount  of  attention.1-25  This  has  been 
driven  by  practical  considerations  and  by  the  severe  difficulties 
inherent  in  experiments  with  real  CWAs.  The  species  thus 
studied  include  Tabun8,22  (GA),  Sarin4-6,10-12’15,17-21’23-25  (GB), 
Soman2,4,10,1 11719,21  (GD),  VX1, 3,7,9,13,14,16  and  their  analogs.  The 
simulant  dimethyl  methylphosphonate  (DMMP)  has  also  been 
studied  computationally.23-31  Of  particular  interest  are  the 
interactions  of  CWAs  with  the  surfaces  of  materials.  This  is  an 
important  issue  in  CWA  detection,  protection  and  remediation 
and  in  understanding  agent  fate,  i.e.,  the  behavior  of  adsorbed 
CWAs  during  exposure  to  the  environment.32 

To  date,  quantum-chemical  work  on  the  interaction  of  CWAs 
with  materials  has  focused  mainly  on  the  adsorption  of  Sarin 
on  clays,  minerals,  and  oxides.19-25  The  interest  here  has  been 
in  determining  the  structure  and  energetics  of  the  adsorbate  and 
in  comparing  observed  and  calculated  infrared  (IR)  vibrational 
spectra.  Much  of  this  work  has  involved  ideal  surfaces  which 
are  free  of  foreign  species.  However  real  oxide  surfaces  exposed 
to  the  environment  are  usually  covered  with  chemisorbed  H20 
(i.e.,  OH  groups)  and  physisorbed  molecular  H20.  This  aspect 
has  received  less  attention.  To  our  knowledge  the  only  ab  initio 
investigations  for  hydroxylated  ionic  metal  oxides  have  been 
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those  of  Michalkova  et  al.: 20,22,26  which  addressed  the  adsorption 
of  Tabun,  Sarin,  and  DMMP  on  MgO  and  CaO  both  with  and 
without  coadsorbed  OH.  It  was  found  that  OH  can  promote  the 
partial  dissociation  of  Sarin  adsorbed  on  MgO  through  the 
elimination  of  an  HF  molecule;  whereas,  the  interaction  of 
Tabun  with  hydroxylated  CaO  is  weak. 

The  present  work  examines  the  behavior,  on  the  hydroxylated 
y-Al203  surface,  of  DMMP,  Sarin,  and  VX  (Figure  1).  There 
are  two  main  objectives.  First  we  wish  to  determine  the  extent 
to  which  hydroxylation  of  the  y-Al203  surface  affects  adsorption, 
in  comparison  to  the  behavior  of  the  OH-free  surface.  Sufficient 
data  are  available  (see  below)  for  adsorption  of  these  reagents 
on  y- AI2O3  to  permit  some  degree  of  comparison  between 
observed  and  calculated  results.  This  effort  will  require  careful 
attention  to  the  details  underlying  the  construction  of  a  suitable 
model  for  the  hydroxylated  y-Al203  surface.  Included  in  this  is 
a  comparison  with  results24  for  adsorption  on  hydroxylated 
amorphous  silica  (a-Si02).  Second,  as  a  step  toward  the 
systematic  understanding  of  issues  relating  to  agent  fate,  we 
are  interested  in  comparing  the  electronic  excitation  of  each 
species  in  the  free  and  adsorbed  states.  This  is  relevant  to 
possible  photochemical  effects  when  CWAs  are  exposed  to 
sunlight.  The  absorption  spectra  of  CWAs  and  simulants  in  the 
near-ultraviolet  (UV)  have  not  been  widely  studied.  Hence,  an 
ability  to  predict  this  property,  and  the  effect  of  adsorption, 
would  be  a  valuable  asset  in  understanding  agent  fate.  This  work 
is  also,  to  our  knowledge,  the  first  computational  study  of  the 
adsorption  of  VX  on  any  material. 
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Figure  1.  Schematic  diagrams  of  the  relevant  molecular  structures 
shown  in  the  optimized  gas-phase  conformations  (see  text).  The  Sarin 
and  VX  enantiomers  were  arbitrarily  chosen. 


After  a  description  of  the  computational  details  and  a  brief 
review  of  the  relevant  experimental  data,  attention  will  focus 
on  the  development  and  testing  of  the  cluster  models  for  the 
hydroxylated  y-Al203  surface.  Next,  a  comparison  of  adsorption 
on  the  hydroxylated  vs  the  OH-free  surface  will  be  presented, 
followed  by  an  examination  of  photoexcitation  in  the  gas-  and 
adsorbed  phases. 

2.  Computational  Methods 

Density  functional  theory  (DFT)  calculations  were  carried 
out  using  the  Amsterdam  density  functional  (ADF  2006  or  2007) 
code.33  An  advantage  of  this  facility  is  the  availability  of 
pseudohydrogen  (PH)  atoms  which  are  hypothetical  neutral 
species  having  a  noninteger  nuclear  charge  and  electron  oc¬ 
cupancy.  As  described  elsewhere,34,35  these  are  useful  for 
terminating  the  coordinatively  unsaturated  atoms  at  the  edges 
of  free-standing  ionic  clusters.  Passivating  the  highly  reactive 
edge  sites  makes  it  possible  to  study  the  interaction  of  large 
and  complex  molecules  such  as  VX  with  small  clusters  without 
incurring  spurious  bond  formation  (see  below).  Pseudo-H 
termination  (PHT)  also  removes  electronic  states,  due  to 
unsaturated  edge  atoms,  in  the  ’’band  gap”  of  the  cluster,  i.e., 
the  gap  between  the  highest  occupied  molecular  orbital  (HOMO) 
and  the  lowest  unoccupied  molecular  orbital  (LUMO).  This 
makes  it  possible  to  treat  problems  related  to  electronic 
excitation  of  adsorbed  species. 

The  basis  sets,  designated  “TZP”,  were  Slater-type  functions 
of  triple-^  quality  with  a  single  polarization  shell  added,  i.e.,  a 
p  shell  for  H  and  a  d  shell  for  all  other  atoms.  These  are  all¬ 
electron  basis  sets  with  no  “frozen  cores”  (cf.  ref  33).  The 
density  functional  consisted  of  the  Vosko— Wilk— Nusair  (VWN) 
form  of  the  local  density  approximation  (LDA)  and  the 
Perdew— Burke— Ernzerhoff  (PBE)  form  of  the  generalized 
gradient  approximation  (GGA)  exchange-correlation  (XC)  po¬ 
tential.  The  PBE  formulation  of  the  GGA— XC  potential  was 
used  because  it  is  compatible  with  analytical  energy  gradients 


in  ADF  and  because  it  has  previously  been  used36-38  in  DFT 
treatments  of  y-Al203  bulk  and  surface  properties. 

For  free  molecules  and  for  adsorption  on  the  small  cluster 
(see  below)  the  convergence  criteria  in  geometry  optimization 
were  1  x  10  6  Hartree  for  total  energy,  1  x  1 0  3  Hartree/A  for 
gradients,  and  1  x  1CT3  A  for  displacements.  The  integration 
precision  parameter  was  set  at  7.0  (tighter  than  the  default  of 
5.0),  and  it  was  verified  that  a  higher  value  (8.0)  had  no 
significant  effect.  For  adsorption  on  the  large  cluster,  which  was 
computationally  demanding,  it  was  necessary  to  relax  the  criteria 
to  the  default  settings  (1  x  1 0  3  Hartree  for  total  energy,  1  x 
1(T2  Hartree/A  for  gradients,  1  x  1CT2  A  for  displacements). 
Further  details  regarding  geometry  optimization  are  given  below. 
The  adsorption  energy  (AZiads)  is  defined  in  the  usual  manner 
as 

A  /-ads =  //(Glister  +  mol)  —  //(cluster)  —  //(mol)  + 

AZZ(BSSE) 

where  the  E' s  are  respectively  the  relaxed  total  energies  of  the 
cluster  with  the  adsorbed  molecule,  of  the  bare  cluster  and  of 
the  free  molecule.  AZi(BSSE)  is  a  counterpoise  correction  for 
basis  set  superposition  error  (BSSE)  which  is  applied  except 
where  noted.  A£ads  <  0  indicates  an  exothermic  reaction. 

Hydrogen  bonding  (H-bonding),  including  C— H-"0  bonds 
between  alkyl  groups  and  O  atoms  in  the  A1203  lattice,  was 
analyzed  using  the  Atoms  in  Molecules  (AIM)  approach39 
developed  by  Bader.40  Critical  points  in  the  total  electron  density 
were  located  using  the  XAIM  code  41  The  vibrational  frequencies 
reported  here  are  the  purely  harmonic  values  with  no  scaling 
correction.  The  appropriate  scaling  factor  in  the  present  case  is 
thought  to  be  close  to  unity.42  The  main  interest  is  in  small  shifts 
in  frequency  caused  by  adsorption  rather  than  in  absolute 
frequencies;  hence,  scaling  should  be  relatively  unimportant. 
Furthermore,  only  the  mid-IR-active  internal  modes  of  the 
adsorbed  species  are  of  interest  here.  Examination  of  the  atomic 
displacements  in  each  of  these  modes  indicates  that,  in  most 
cases,  there  is  little  or  no  contribution  from  motion  of  atoms  in 
the  A1203  cluster  itself.  Hence,  the  frequency  results  are  reliable 
even  though  only  partial  optimization  of  the  cluster  geometry 
is  performed  (see  below). 

It  is  known43  that  proper  description  of  excited  electronic 
states  requires  the  inclusion  of  diffuse  functions  in  the  basis 
sets.  In  the  present  case,  QZ3P+1  basis  sets  and  time-dependent 
density  functional  theory  (TDDFT)  were  employed  when 
treating  electronic  excitation.  These  basis  sets  are  of  quadruple-^ 
quality  with  three  polarization  and  one  diffuse  function  added. 
The  TDDFT  studies  were  done  as  single-point  calculations  for 
stmctures  optimized  using  time-independent  DFT  and  TZP  basis 
sets.  The  XC  potential  used  in  TDDFT  was  based  on  the 
statistical  averaging  of  orbital  potentials  (SAOP)  method44  which 
has  the  correct  —  1/r  dependence  for  r  — *  Tests  were  done  to 
verify  that  the  computed  excitation  energies  are  independent 
of  the  technical  aspects  of  the  TDDFT  calculation  (e.g., 
integration  accuracy,  SCF  convergence  criteria,  etc.).  For 
studying  the  near-UV  spectra  of  adsorbed  species,  QZ3P+1 
basis  sets  were  used  for  the  molecule  and  for  the  Al(— O— )3 
site  at  which  adsorption  occurs  (see  below).  In  order  to  reduce 
the  computational  time,  TZP  basis  sets  were  used  for  the  rest 
of  the  cluster.  No  significant  difference  was  seen  when  the 
Al(— O— )3  adsorption  site  was  also  treated  with  TZP  basis  sets. 
Some  TDDFT  calculations  were  also  done  for  free  molecules 
using  the  Gaussian  03  code45  with  the  B3LYP  functional  and 
6-311++G(d,p)  basis  sets.  These  used  “Tight”  convergence 
criteria  in  geometry  optimizations  (with  time-independent  DFT) 
and  an  “UltraFine”  integration  grid. 


Adsorption  of  DMMP,  Sarin,  and  VX 

3.  Experimental  Background 

This  section  briefly  summarizes  the  relevant  experimental  data 
for  DMMP,  Sarin,  and  VX  both  free  and  adsorbed  on  y-ALOi 
and  for  the  photochemical  properties  in  the  near-UV. 

3.1.  Free  Molecules.  The  conformations  of  gas-phase  DMMP 
and  Sarin  have  been  studied  both  experimentally12,27  (using 
microwave  spectroscopy)  and  theoretically .4.10, 12,27,46,47  Gas-phase 
vibrational  spectra  (which  differ  somewhat  from  those  for  the 
liquid)  have  been  reported  for  DMMP  (refs  48—50)  and  Sarin 
(ref  51);  see  also  ref  23.  During  geometry  optimization,  here  as 
before, 23-25  DMMP  and  Sarin  relaxed  into  their  lowest-energy 
conformations  (Figure  1)  which  concur  with  those  reported 
previously.10,27  For  Sarin,  calculations  of  the  relaxed  total  energy 
vs  the  angle  of  rotation  of  the  (CH3)2(H)C—  group  about  the 
C-OP  bond  were  done  in  10°  steps.  The  results  (not  shown) 
were  very  similar  to  those  obtained  previously10  at  the  B3LYP/ 
6-31G(d,p)  level.  In  Sarin  and  VX,  the  P  atom  is  an  optically 
active  center.  In  either  case,  only  one  arbitrarily  chosen 
enantiomer  (that  shown  in  Figure  1)  was  considered  here. 

Rotational  isomerization  in  VX  and  related  species  has  been 
studied  in  dilute  solutions  using  IR  spectroscopy52,53  and 
semiempirical  theory.53  The  only  isomerization  detected  involves 
hindered  rotation  about  the  CH3CH2O— P  bond.  However,  VX 
in  solution  assumes  a  more  compact  geometry  than  in  the  vapor 
phase,  for  which  other  forms  of  isomerization  may  be  possible. 
Rotational  isomerization  in  isolated  VX  is  discussed  in  more 
detail  below  and  in  the  Supporting  Information  in  connection 
with  steric  effects  in  adsorption.  An  ab  initio  computational 
study1  for  vapor-phase  VX  led  to  the  optimized  structure  shown 
in  Figure  lc.  A  similar  study,  at  the  B3LYP/6-31+G(d)  level, 
has  been  reported7  for  an  analog  of  VX  (termed  ”VX-Me”)  in 
which  the  ethoxy  group  is  replaced  by  a  methoxy  and  the 
isopropyl  groups  by  methyl  groups.  Although  different  enan¬ 
tiomers  were  considered  in  the  two  studies,  comparison  of  the 
torsion  (dihedral)  angles  shows  that  the  optimized  structures 
are  similar.  Vibrational  spectra  are  available54-56  for  VX  with, 
in  the  case  of  Raman  data,56  detailed  mode  assignments. 

3.2.  Adsorption  on  v-ALOv  The  interaction  of  DMMP  with 
high-surface-area  (HSA)  y-AFOs  powder  has  been  studied  by 
Mitchell  et  al.57  Results  have  also  been  reported  for  adsorption 
on  an  unspecified  form  of  AI2O3  HSA  powder,58  on  oxide  layers 
on  polycrystalline  Al  films59  (believed  to  be  y-Al203),  and  on 
an  oxide  layer  on  an  Al(lll)  surface.60  For  OH-free  surfaces 
at  room  temperature,  anondissociative  dative-bonding  interaction57, 59 
occurs  between  the  coordinatively  unsaturated  tetrahedral  All!)/) 
Lewis  acid  site  and  the  O  atom  of  the  P=0  group.  Surface  OH, 
if  present,  may  also  be  involved  through  H-bonding  to  the  P=0 
group  and/or  to  an  —  OCH3  group.  Higher  temperature  promotes 
a  nucleophilic  attack  by  the  O  atom  of  an  OH  site  on  the  P 
atom,  causing  the  release  of  CH3OH  and  the  formation  of  a 
bridging  structure.  For  nondissociative  adsorption572,592  only  the 
v(P=0)  mode,  at  1216  cm-1,  exhibits  a  substantial  shift  from 
the  gas-phase  value48-50  (1276  cm-1).  Other  modes  appear  at 
energies  close  to  the  gas-phase  values.  Reaction  with  OH  affects 
modes  related  to  the  CH30—  groups  and  further  increases  the 
red-shift  in  v(P=0).  At  still  higher  temperatures59  loss  of  the 
second  CH30  occurs  with  formation  of  an  (— Al— O— )3PCH3 
surface  species.  None  of  the  data  show  any  clear  indication  of 
shifts  in  the  v(0— H)  modes  of  surface  hydroxyls  due  to 
H-bonding  to  DMMP.  However,  such  bonding  is  known61  to 
result  in  a  very  broad  v(0— H)  band  which  can  be  difficult  to 
detect. 

The  adsorption  of  Sarin  on  HSA  y-Al203  powder  has  been 
studied  experimentally  using  IR  spectroscopy.62,63  Initial  adsorp- 
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tion  occurs  via  a  strong  Al •  0=P  interaction  which,  on  an 
hydroxylated  surface,  is  followed  by  H2O  elimination,  transfer 
of  F  to  an  adjacent  Al  site,  and  formation  of  a  bridging  structure 
analogous  to  that  described  above  for  DMMP.  For  nondisso¬ 
ciative  adsorption  the  IR  data  show,  as  the  only  major  effect,  a 
shift  of  Av(P=0)  =  —63  cm-1  relative  to  the  gas-phase  value51 
of  1308  cm-1.  As  in  the  case  of  DMMP,  subsequent  reaction 
with  OH  groups  leads  to  more  substantial  changes  in  the  IR 
spectra. 

Experimental  data  for  VX  adsorption  on  A1203  are  provided 
in  the  NMR  studies  of  Wagner  et  al. 64,65  which  involve  the  use 
of  both  aerogel-prepared  A1203  nanoparticles  and  y-Al203  HSA 
powders.  For  high  VX  loading,  complete  reaction  of  the  entire 
nanoparticle  occurs  to  form  an  aluminophosphonate  complex. 
There  is  also  evidence  that  H20  and/or  OH  on  the  HSA  powder 
promotes  reaction  with  VX  leading  to  a  bidentate  phosphonate 
structure. 

3.3.  UV  Spectra  and  Photochemical  Properties.  A  sub¬ 
stantial  amount  of  work  has  been  done  on  the  photocatalytic 
destruction  of  DMMP  using  TKL  and  other  metal  oxides  (see, 
e.g.,  refs  66  and  67  and  works  cited).  These  photoeffects  are 
initiated  by  the  excitation  of  electron— hole  pairs  in  the  oxide. 
In  the  present  work,  we  are  concerned  with  the  photoexcitation 
of  CWAs  themselves  when  adsorbed  on  a  photochemically  inert 
oxide.  With  a  bulk  band  gap  of  ~7.2  eV  (ref  68),  y-Al203  is 
transparent  throughout  the  photochemically  active  part  of  the 
terrestrial  solar  spectrum69  (TSS),  the  intensity  of  which  is 
negligible  above  about  4.5  eV  (see  below). 

The  UV  spectra  of  DMMP  and  Sarin,  in  the  vapor  phase  or 
dissolved  in  cyclohexane  or  ethanol,  have  been  reported.70-72 
Both  are  weak  absorbers  at  wavelengths  longer  than  ~200  nm 
(energies  below  ~6.2  eV).  In  cyclohexane,70  DMMP  shows  a 
weak  band  (possibly  due  to  an  impurity)  at  205  nm  (6.1  eV), 
and  in  the  vapor  phase71  there  appears  to  be  a  stronger  band  (or 
shoulder)  at  ~160  nm  (7.8  eV).  This  shoulder  lies  on  the  edge 
of  what  appears  to  be  a  strong  absorption  continuum.  Sarin  in 
the  vapor  phase  shows  no  absorption  for  X  >  200  nm  but  in 
ethanol  exhibits  a  weak  band  at  about  204  nm  which  may  be 
due  to  an  impurity.72  Aschmann  et  al.73  observed  no  significant 
photochemical  decomposition  of  gas-phase  DMMP  during 
irradiation  at  X  >  300  nm  ( hv  <4.1  eV).  Rauk  et  al.15  have 
calculated  the  lower-lying  electronic  transitions  of  a  Sarin  analog 
with  CH3  replacing  (CH3)2CH.  At  the  MP2/6-31+G(d)  level  the 
excitation  threshold  is  found  at  about  9.9  eV,  somewhat  higher 
than  that  computed  here  for  Sarin  (see  below). 

VX  dissolved  in  cyclohexane  was  found70  to  be  a  somewhat 
stronger  UV  absorber  than  DMMP  or  Sarin;  although,  no  actual 
absorption  peaks  were  seen  for  X  >  200  nm.  Zuo  et  al.74  studied 
VX  dissolved  in  ethanol  and  found  absorption  maxima  at  209, 
238,  and  266  nm  (5.93,  5.21,  and  4.66  eV).  These  were  analyzed 
(for  the  gas-phase  molecule)  at  the  restricted  Hartree— Fock 
(RHF)  level  using  the  configuration  interaction  singles43  ap¬ 
proach  with  STO-3G  basis  sets.  Irradiation  of  VX  droplets  in 
air  at  253.7  nm  (4.89  eV)  led  to  a  complex  series  of  molecular 
rearrangements  and  also  to  reaction  with  photogenerated  ozone. 
The  reported  UV  spectra70,74  for  VX  differ  qualitatively; 
although,  the  effect  of  the  solvent  (cyclohexane  vs  ethanol)  is 
uncertain. 

4.  Construction  of  y-ALO,  Surface  Models 

This  section  describes  the  construction  and  evaluation  of  the 
cluster  models  for  the  hydroxylated  surface.  It  will  be  shown 
that  the  models  give  a  good  description  of  both  the  chemical 
and  the  electronic  structure  of  the  hydroxylated  y-Al203  surface 
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Figure  2.  Models  for  (a)  the  defective-spinel  and  (b)  the  ideal  spinel 
y-Al2C)3  lattices.  Blue,  green,  and  red  spheres  are  respectively  Al(7’</), 
A1(0;,),  and  O.  In  (a)  the  small  black  spheres  represent  A1(0/,)  vacancies. 
Different  surface  planes  are  identified. 


Figure  3.  Models  (not  to  scale)  for  the  (a)  AlgOi2  and  (b)  AI32O48 
clusters.  For  clarity,  the  pseudo-H  atoms  (H1/2,  etc.)  are  shown  only 
for  (a).  Examples  of  the  different  O-atom  types  (Oa  and  Ob)  are  labeled. 
At  the  “surface”,  an  Oa  is  bonded  to  two  A1(0/,)  and  one  Al(7)j),  and 
an  Ob  is  bonded  to  three  A1(0/,).  The  “1”  labels  the  central  Al(Trf)  site, 
and  “2”  and  “4”  label  atoms  according  to  their  distance  from  this  site. 
Site  “3”  would  be  a  subsurface  A1(0/,).  Thus  placing  an  OH  on  the 
central  A1  and  H  on  the  nearest  Ob  gives  a  1,4  configuration. 


and  that  there  is  no  pronounced  dependence  of  the  properties 
of  the  OH  groups  on  the  size  of  the  PH-terminated  clusters. 
These  results  are  general  and  will  be  useful  in  other  studies  of 
the  hydroxylated  y-Al203  surface. 

4.1.  OH-Free  Surface.  The  bulk  and  surface  structures  of 
y-AhOg  are  controversial  and  are  discussed  briefly  in  refs  23 
and  25  together  with  a  description  of  cluster  models  for  the 
OH-free  surface.  The  controversy  centers  on  whether  a  defec¬ 
tive-spinel  or  a  nonspinel  model  is  more  appropriate,  and  a 
recent  study75  (which  reviews  earlier  work)  suggests  that 
y-Al203  should  actually  be  seen  as  an  amorphous  random 
network  of  Al— O  bonds.  The  model  for  the  bulk  structure 
determines  the  structure  and  stability  of  the  different  surface 
planes  which,  in  turn,  influences  the  identity  and  the  distribution 
of  reactive  surface  sites  (both  OH-free  and  hydroxylated). 

In  the  present  work,  as  before,  the  point  of  departure  is  the 
defective-spinel  model  for  bulk  y-Al203  described  by  Pinto  et 
al. 76,77  and  shown78  in  Figure  2a.  Further  details  are  given  in 
the  Supporting  Information.  Following  convention,  the  y-Al203 
surface  planes  are  labeled  according  to  the  corresponding  planes 
of  the  spinel  lattice  from  which  they  are  derived.  Thus  the  (010) 
plane  of  the  monoclinic  defective-spinel  unit  cell  in  Figure  2a 
corresponds  to  the  (111)  plane  of  the  cubic  spinel  lattice,  Figure 
2b.  The  most  stable  surface  in  this  model  (the  (lll)a  in  the 
spinel  labeling)  is  formed  by  ’’cleaving”  on  the  layer  of 
vacancies  shown  in  Figure  2a.  This  surface  consists  of  singly 
unsaturated  tetrahedral  Al (Tj)  sites  and  two  types  of  O  sites. 
One  site  (termed  ”Oa”  in  the  present  work)  is  3-fold  coordinated 
in  the  bulk,  as  a  result  of  being  adjacent  to  a  vacant  octahedral 


A1(0/,)  site,  and  remains  so  at  the  surface.  An  Oa  has  one  AKTj) 
and  two  A1(0/,)  nearest-neighbors.  The  other  O  site  (termed 
“Ob”)  is  4-fold  coordinated  in  the  bulk  but  is  missing  its  Al(2j/) 
nearest-neighbor  at  the  surface  and  thus  has  three  A1(0;,)  nearest- 
neighbors.  There  are  no  unsaturated  Al {Oh)  sites  on  the  (111 )a 
surface. 

Figure  3  shows  two  clusters  cut  from  the  (lll)a  surface 
shown  in  Figure  2a.  These  are  of  stoichiometry  A180i2  and 
Al32048  and  are  used  with  PHT  as  described  in  the  Supporting 
Information.  In  both  cases,  the  reactive  cation  surface  site,  in 
the  absence  of  OH,  is  an  Al(7^)  which  is  known  from 
experiment57,59,62  to  be  important  in  the  adsorption  of  DMMP 
and  Sarin.  In  the  first  step  in  geometry  optimization  for  OH- 
free  surfaces,  the  Al  and  O  atoms  in  the  cluster  were  frozen  in 
the  ideal  bulk-lattice  positions  while  the  PHs  were  all  allowed 
to  vary.  In  subsequent  optimizations,  most  of  the  PH,  Al,  and 
O  atoms  remained  fixed  (see  below). 

4.2.  Hydroxylated  Surface.  4.2.1.  Chemical  Structure.  The 
structure  of  the  hydroxylated  y-Al203  surface  is  complex  and 
not  yet  fully  understood.  Many  chemically  distinct  OH  sites 
are  possible,  depending  on  the  surface  structure  and  treatment. 
Empirical  models  have  been  proposed,79-83  based  on  analyses 
of  the  IR  spectra  of  HS  A  powders  and  the  response  to  chemical 
treatment,  which  correlate  the  v(0— H)  stretching  frequency  with 
the  structure  of  the  OH  site.  Recently  these  models  have  been 
further  investigated  using  quantum-chemical  methods.84-86  The 
general  consensus  is  that  v(0~ H)  falls  in  approximately  the 
3760—3800  cm-1  range  for  type  I  sites  (OH  bonded  to  a  single 
surface  Al),  3700—3750  cm-1  for  type  II  sites  (OH  bridging 
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two  Al’s),  and  3635—3700  cm-1  for  type  III  sites  (OH  bonded 
to  three  Al‘s).  For  each  type  of  site,  various  combinations  of 
AKO/J  and  A I  (7),)  are  possible.  For  example,  a  type  II  OH 
bridging  two  A1  (Oh)  sites  is  chemically  distinct  from  one 
bridging  an  A1  (Oh)  and  an  Al( Td). 

The  (111  )a  surface  described  above  differs  from  both  the 
( 1 1 1 ) A  and  (111  )B  surfaces  typically  invoked  in  discussions79,82 
of  the  spinel-like  y-Al203  lattice.  Figure  2b  shows  the  ideal- 
spinel  y-AFOi  lattice  which  involves  no  A1  vacancies  and  is 
therefore  nonstoichiometric  (see  Supporting  Information).  Cleav¬ 
ing  on  the  plane  indicated  generates  the  (lll)A  and  (111 )B 
surfaces.  When  fully  reacted  with  H20,  the  former  exhibits  type 
I  AK7)/)— OH  and  type  II  A\(Ok)~0(\\)-A\(Ttl)  sites  in  a  1:3 
ratio.  The  latter  involves  type  II  A1(0/,)— O(H)— A1(0/,)  and  type 
III  [Al(0/,)]3— OH  sites  in  a  3:1  ratio.  These  points  are  discussed 
more  fully  in  refs  79  and  82.  When  the  optimal  distribution  of 
A1(0/,)  vacancies  is  considered  in  the  defective-spinel  model,76,77 
the  (111  )cz  plane  shown  in  Figure  2a  becomes  the  lowest-energy 
surface  since  it  cuts  through  a  high  density  of  vacancies.  On 
this  surface,  a  type  I  AllT,/)— OH  site  and  type-III 
[AKO^MATOJ-OaH  or  [Al(0/,)]3— ObH  sites  are  formed. 
Henceforth,  these  type  III  sites  will  be  termed  “type  Ilia”  and 
“type  Illb.  Two-dimensionally  periodic  slab  (2-DPS)  calcula¬ 
tions77  for  the  dissociative  adsorption  of  H20  show  little 
dependence  of  AE^  (at  most  ~5%)  on  which  O  site  is  chosen 
for  the  H  atom  (see  below). 

In  order  to  construct  a  viable  computational  model  for 
adsorption  on  the  hydroxylated  surface,  one  which  can  be 
compared  with  data  for  HSA  powders,  it  is  necessary  to 
determine  which  OH  sites  are  important.  For  DMMP,  Mitchell 
et  al.57a  found  that  the  OH  sites  removed  in  the  initial  adsorption 
at  room  temperature  correspond  to  a  loss  of  intensity  in  the 
range  of  ~3700— 3750  cm-1.  Similar  data  have  been  reported59 
for  DMMP  adsorption  on  A1203  films  on  polycrystalline  Al. 
These  were  thought  to  be  of  the  y  form  and  free  of  physisorbed 
H20.  Here  adsorption  led  to  a  loss  of  intensity  in  the  3630—3650 
cuT1  range.  For  Sarin  the  reactive  OH  group  corresponds62  to 
a  v(0~ H)  of  3700  cm-1.  These  results  all  exclude  type  I  sites 
as  the  active  species,  but  type  II  and  type  III  remain  as 
possibilities.  Type  I  sites  may  be  present  as  spectators  but,  in 
view  of  the  IR  data,  they  do  not  participate  in  the  initial 
adsorption  of  DMMP  or  Sarin. 

The  models  in  Figure  3  were  hydroxylated  by  dissociatively 
adsorbing  H20  (see  below).  An  OH  was  placed  at  the  center  of 
the  surface,  to  form  a  type  I  Al(Td)— OH,  and  an  H  at  either  an 
Oa  or  Ob  site  to  form  a  type  Ilia  or  Illb  OH  site.  The  discussion 
thus  far  has  been  entirely  in  terms  of  the  defective-spinel  model 
for  y- A1203.  Examining  the  nonspinel  model84  one  finds  that 
type  Illb  sites  also  form  on  the  (100)  and  (111)  surfaces  and 
that  this  is  the  only  type  III  site  (termed  “HO-/r3-AlVI”  in  ref 
84)  on  the  low-index  surfaces  of  this  model.  The  type  I 
Al(7rf)— OH  (termed  “HO-z^-A^v”  in  ref  84)  is  also  found  on 
the  (110)  nonspinel  surface.  Hence  the  results  obtained  here 
should  be  generally  applicable  to  both  y-Al203  models.  The 
OH  sites  treated  here  consist  of  the  basic  type  I  and  the  acidic 
type  III,  and  it  is  found  that  OH" •  0=P  bonding  clearly  favors 
the  latter  (see  below).  This  limiting  behavior  then  provides  a 
firm  basis  for  estimating  (in  lieu  of  actual  results)  the  behavior 
of  intermediate  sites  such  as  the  type  II  found  on  other  surfaces 
of  both  the  spinel  and  nonspinel  models.  Likewise,  the  present 
models  consider  only  the  Al (Td)  site  on  the  OH-free  surface. 
The  nonspinel  model846  also  involves  unsaturated  Al  (Oh)  surface 


sites,  but  these  are  found  to  be  lower  in  Lewis-acidity  and  are 
thus  expected25  to  be  less  active  in  dative  bonding  to  the  P=0 
group. 

Geometry  optimization  was  performed  for  the  hydroxylated 
models.  These  were  partial  optimizations  which  included  all 
Al  atoms  in  the  ’’surface”  (one  for  A180i2  and  seven  for  Al32048), 
the  three  O  nearest-neighbors  of  the  type  I  A\(Td)— OH  and  the 
Oa  or  Ob  atom  to  which  the  H  was  added.  The  adsorbed  H  and 
OH  were  also  included  in  the  optimization.  The  rest  of  the  Al 
and  O  atoms  remained  fixed  in  the  bulk-lattice  positions,  and 
the  PHs  were  frozen  in  the  previously  determined  optimized 
positions. 

4.2.2.  Formation  via  Dissociative  Adsorption  ofH20.  This 
section  examines  A£ads  for  the  chemisorption  of  H20  on  the 
OH-free  surface  which  is  a  useful  quantitative  test.  The 
energetics  of  the  process  whereby  the  OH-free  cluster  is 
converted  to  the  proposed  model  hydroxylated  surface  should 
be  consistent  with  results  for  the  chemisorption  of  HzO.  The 
related  issue  of  the  physisorption  of  molecular  H20,  in 
comparison  to  2-DPS  results,  is  discussed  elsewhere25  as  a 
quantitative  test  of  the  OH-free  model  surface. 

Previous  2-DPS  results77  for  dissociative  adsorption  of  H20 
on  the  (1 1  l)a  surface  showed  little  dependence  of  A Elis  on  the 
placement  of  the  H  atom,  relative  to  the  type  I  A\(Td)  site,  or 
on  whether  a  type  Ilia  or  Illb  site  was  formed.  A  A Z?ads  of  —36.2 
to  —38.0  kcal/mol  was  found  for  all  of  the  configurations 
considered.  In  the  present  work,  2-DPS  calculations25  were  done 
for  dissociative  adsorption  in  a  1,4  configuration  with  H  at  a 
type  Illb  site  and  in  a  1,2  configuration  with  H  at  a  type  Ilia. 
After  correction  for  BSSE,  the  former  gave  AEads  =  —33.9  kcal/ 
mol,  and  the  latter  gave  —34.9  kcal/mol.  In  this  notation  “1” 
refers  to  the  Al  site  with  the  type  I  OH  group,  and  “2”  or  “4” 
refers  to  placement  of  the  H  atom  on  a  first-  or  third-nearest- 
neighbor  O  atom.  The  second-nearest-neighbor  of  an  A\(Td)  is 
an  Al(Oft);  hence,  no  1,3  OH  configuration  is  formed.  These 
results  will  now  be  compared  with  those  for  the  clusters. 

For  the  AlgOi2  cluster  (Figure  3a),  an  OH  is  placed  on  the 
Al/Tj)  site  and  an  H  on  the  Ob  atom  in  a  1,4  configuration.  The 
latter  is  a  type  Illb  site  if  the  H1/2  PH  is  counted  as  the  third 
Al  (Oh)  neighbor.  This  is  justified  in  that  the  properties  of  this 
site  (see  below)  compare  well  with  those  of  a  1,4  type  Illb  site 
on  the  larger  Al32048  cluster.  This  structure  gave  AFads  =  —33.8 
kcal/mol  which  compares  well  with  the  present  2-DPS  result 
(—33.9  kcal/mol).  Hence,  even  the  small  A180i2  cluster,  with 
PHT,  can  give  good  results  for  the  chemisorption  of  H20. 
However,  placing  the  H  on  the  Oa  site  gives  a  1,2  configuration 
with  AEais  =  —29.9  kcal/mol  which  differs  somewhat  from  the 
present  2-DPS  result  (—34.9  kcal/mol). 

Similar  calculations  were  done  for  the  Al32048  cluster  (Figure 
3b).  The  1,4  configuration  with  H  on  an  Ob  site  gave  Aitads  = 
—37.3  kcal/mol  which  is  fairly  close  to  the  2-DPS  result  (—33.9 
kcal/mol).  The  optimized  type  III  Ob— Hads  distance  was  0.978 
A.  For  the  type  I  Al— OH,  the  Al— O  and  O— H  distances  were 
1.758  and  0.971  A,  respectively,  and  the  Al— O— H  bond  angle 
was  112°.  The  larger  Al32048  cluster  allows  more  flexibility  in 
the  choice  of  the  site  for  the  H  atom.  Another  1,4  configuration, 
but  one  involving  an  Oa  site,  gave  A/s^s  =  —28.9  kcal/mol 
which  again  differs  somewhat  from  the  2-DPS  result.  This  result 
is,  however,  close  to  the  value  (—29.9  kcal/mol)  obtained  above 
for  a  1,2  type  Ilia  configuration  on  the  A180i2  cluster. 

To  summarize  the  cluster  results,  the  1,2  and  1,4  type  Ilia 
structures  give  very  similar  energies,  whereas  1,4  type  Ilia  and 
1 ,4  type  Illb  configurations  give  somewhat  different  values  with 
the  latter  being  in  better  agreement  with  2-DPS  results.  Thus, 
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Figure  4.  (a)  Computed  DOS  for  the  hydroxylated  AI32O48  cluster 
shown  in  Figure  3b.  The  zero  of  energy  is  the  VBM.  A  bonding  state 
derived  from  the  type-IIIb  OH  site  is  indicated.  The  DOS  has  been 
broadened  with  an  arbitrarily  chosen  Gaussian  width  of  0.27  eV.  (b) 
SXPS  data  (ref  87  with  permission)  for  an  hydroxylated  AI2O3  film  on 
A1  (110).  The  energy  scale  has  been  shifted  by  1.0  eV  from  that  given 
in  the  original  reference  to  place  the  VBM  at  the  foot  of  the  strong 
peak  so  as  to  be  consistent  with  (a).  The  feature  assigned  in  ref  87  to 
OH  is  marked.  “O  2s”  and  "O  2p”  indicate  the  dominant  orbital 
composition  of  the  two  SXPS  bands. 

Aliads  is  nearly  independent  of  the  distance  between  the  type  I 
and  III  sites,  in  agreement  with  the  2-DPS  models,  but  there  is 
a  dependence  on  whether  a  type  Ilia  or  Illb  site  is  formed  which 
is  not  seen  for  the  2-DPS. 

4.2.3.  Surface  Electronic  Structure.  The  discussion  thus  far 
has  focused  on  the  accuracy  of  the  cluster  model  with  respect 
to  surface-chemical  properties.  Another  test  is  based  on  the 
electronic  structure.  In  the  following,  reference  will  be  made 
to  the  “valence  band”  (VB)  of  the  cluster  and  to  the  “valence 
band  maximum”  (VBM),  neglecting  the  fact  that  the  cluster 
orbitals  constitute  a  finite  (albeit  closely  spaced)  set  of  discrete 
states.  Soft-X-ray  photoemission  spectroscopy  (SXPS)  data87 
for  hydroxylated  AI2O3  films  on  an  A1  (110)  surface  (Figure 
4b)  show  a  bonding  state,  for  an  unknown  type  of  OH  site,  at 
about  9  eV  below  the  VBM.  Figure  4a  shows  the  density  of 
states  (DOS)  computed  for  the  1,4  type  Illb  configuration  of 
the  AI32O48  cluster  (Figure  3b).  A  deep-lying  state,  due  to  the 
type  Illb  OH  site,  is  seen  at  10.6  eV  below  the  VBM.  Results 
(not  shown)  for  the  1 ,4  type  Ilia  configuration  are  similar,  with 
the  type  Ilia  OH  feature  at  8.8  eV  below  the  VBM.  For  the 
1 ,4-hydroxylated  Al80i2  cluster  the  type  Illb  OH  peak  falls  at 
9.1  eV  below  the  VBM. 

To  investigate  this  issue  further,  DOS  results  (not  shown) 
were  obtained  for  the  present  2-DPS  calculations.  These  gave 
a  type  Illb  OH  peak  at  9.6  eV  below  the  VBM  in  good 
agreement  with  experiment  and  with  the  cluster  calculations. 
The  type  I  OH  gave  a  broad  peak,  centered  at  ~5.8  eV  below 
the  VBM,  due  to  the  AlO— H  bond  and  a  narrower  peak,  just 


TABLE  1:  Deprotonation  Energies  for  the  Hydroxylated 
Clusters" 

OH  type  A18012  AI32O48 

I  (-1600)'' 

Ilia  1228  1161 

Illb  1244  1199 

a  Energies  are  given  in  kJ/mol  to  conform  with  refs  88  and  89 
and  have  been  corrected  for  a  vibrational  zero-point  energy  of  30 
kJ/mol  associated  with  the  O— H  stretching  and  bending  modes  (see 
text).  All  values  are  positive  (endothermic).  The  smaller  the  A.Edp 
the  more  acidic  the  OH  group.  Both  clusters  are  pseudo-H 

terminated.  b  The  deprotonated  type  I  site  is  metastable  (see  text). 

This  result  is  an  estimate. 

below  the  VBM,  due  to  the  Al— OH  bond.  For  the  AI32O48 
(Al80i2)  cluster  the  corresponding  type  I  states  lie  at  ~5.9  (5.4) 
and  1.2  (1.0)  eV  below  the  VBM.  The  type  I  OH  features  are 
not  readily  observable  in  SXPS  due  to  the  strong  emission  from 
the  AI3O3  VB;  however,  both  clusters  give  DOS  results  in  good 
agreement  with  the  2-DPS  and  with  the  available  experimental 
data. 

4.2.4.  Br0nsted  Acidity.  The  acidity  of  the  different  OH  types 
is  important  in  reactivity.  Nortier  et  al.88  have  estimated  the 
Brpnsted  acidity  of  OH  on  d-  and  O-AI2O3  by  computing  the 
deprotonation  energy  (A£dp)  at  the  RHF  level  using  small  cluster 
models  and  effective-core  pseudopotentials.  For  deprotonation 
of  a  species  XOH,  AEdp  is  defined  as  F(XO  )  +  E(H+)  — 
E(XOH)  where  XO  and  H+  are  at  infinite  separation  and 
experience  no  Coulombic  attraction.  Aiidp  is  positive  (endo¬ 
thermic)  and  decreases  with  increasing  acidity.  For  AI2O3,  A£dp 
decreases  in  the  order  type  I  >  type  II  >  type  III.  The  type  Ilia 
site  is  found88  to  be  the  most  acidic,  more  so  than  the  type  Illb. 
The  type  I  N\(Td)— OH  exhibits  the  second-lowest  acidity,  with 
only  the  type  I  A1(0/,)— OH  site  being  the  less  acidic.  A  similar 
ordering  of  acidities  was  also  obtained  in  another  study89  of 
small  clusters  at  the  RHF/3-21G  level. 

Deprotonation  energies  (Table  1)  were  obtained  here  using 
the  DFT  methods  and  hydroxylated  clusters  described  above. 
These  have  not  been  corrected  for  BSSE,  but  based  on  results 
obtained  above,  including  A E  (BSSE)  would  reduce  all  Aftdp 
values  by  only  ~5  kJ/mol.  A  correction  has  been  applied  for 
the  zero-point  energies  of  the  O— H  stretching  and  bending 
modes  which  are  lost  upon  deprotonation.  This  reduces  the 
purely  electronic  A£dp  by  ~30  kJ/mol.  The  ordering  of  acidities 
in  Table  1  agrees  with  that  in  previous  work;  although,  the 
difference  between  type  Ilia  and  Illb  sites  is  not  as  large.  Also, 
the  AlgOi2  results,  although  larger  than  the  corresponding 
AI32O48  values,  are  reasonably  close  to  those  for  the  larger 
cluster. 

The  calculations  differ  somewhat  from  those  of  Nortier  et 
al.88  in  that  two  types  of  OH  sites  are  present  when  one  or  the 
other  is  deprotonated.  For  example,  when  deprotonating  a  type 
I  OH  there  is  also  a  nearby  type  III  OH  which  retains  its  proton. 
This  situation  is  more  relevant  to  the  present  investigation  than 
is  an  isolated  OH.  It  was  found  here  that  the  deprotonated  type 
I  site  is  metastable  on  either  cluster.  Relaxation  leads,  ultimately, 
to  migration  of  the  H  atom  from  the  type  Illb  site  to  the 
AlfT,/)— O  to  reform  the  type  I  OH  at  the  expense  of  the  type 
Illb.  This  is  a  clear  indication  of  the  lower  type  I  OH  acidity 
since  the  deprotonated  type  I  site  lies  higher  in  energy  than  the 
deprotonated  type  III.  An  estimate  of  AEdp  ~  1600  kJ/mol  for 
the  type  I  OH  is  obtained  from  the  energy  of  the  metastable 
configuration  formed  after  the  initial  large  relaxation  but  before 
the  start  of  pronounced  H-atom  migration. 
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TABLE  2:  Structural  and  Energetic  Results  for  Nondissociative  Adsorption  of  DMMP,  Sarin,  and  VX  on  PH-Terminated 
Clusters0 


cluster 

adsorbate 

Aiiads 

/-(Al— OP) 

r(A10=P)  b 

Z(A1— 0=P) 

Av(P  =  O) 

A18012  (OH-free) 

DMMP 

-53.0 

1.819 

1.523  (1.490) 

149 

—62c 

Sarin 

-51.5 

1.838 

1.521  (1.481) 

132 

-103d 

VX 

-51.6 

1.814 

1.537  (1.497) 

137 

-95e 

A132048  (OH-free) 

DMMP 

-41.4 

1.860 

1.516 

148 

Sarin 

-39.8 

1.863 

1.513 

130 

VX 

-39.3 

1.860 

1.534 

136 

"In  all  cases,  adsorption  is  via  AI(7),)"-0=P  dative  bond  formation.  Energies  are  in  kcal/mol,  bond  lengths  in  Angstroms,  bond  angles  in 
degrees  and  frequencies  in  cnT1.  The  Av(P=0)  values  are  based  on  the  harmonic  values  with  no  scaling  correction.  Av(P=0)  is  v(P=0) 
calculated  for  the  adsorbed  species  minus  that  calculated  for  the  gas-phase.  The  DMMP  and  Sarin  results  repeat  those  in  ref  25.  All  energies 
are  BSSE-corrected.  b  Numbers  in  parentheses  are  the  P=0  bond  lengths  computed  for  the  free  molecules.  c  The  experimental  result  is  —60 
cm-1  (ref  57),  based  on  a  gas-phase  v(P=0)  of  1276  cm-1  (refs  48—50).  For  adsorbed  DMMP  the  computed  P=0  stretch  is  strongly  coupled 
to  modes  at  1161  and  1174  cnT1,  and  v(P=0)  is  taken  as  the  average.  rfThe  experimental  result  is  —63  cm-1  (ref  62),  based  on  a  gas-phase 
v(P=0)  of  1308  cnT1  (ref  51).  'No  corresponding  experimental  data  are  available. 


5.  Results  and  Discussion 

5.1.  Adsorption  on  OH-Free  Surfaces.  The  adsorption 
energies  and  geometries  for  DMMP  and  Sarin  on  OH-free 
y-Al203,  computed  using  PH-terminated  clusters,  are  given 
elsewhere25  and  are  repeated  in  Table  2  together  with  new  results 
for  VX.  The  dependence  of  A Eads  on  cluster  size  and  on 
adsorbate  has  been  discussed25  in  terms  of  the  energy  difference 
between  the  donor  level  of  the  free  molecule  and  the  acceptor 
level  associated  with  surface  Al(7)/)  site.  For  the  defective-spinel 
model76-77  for  y-Al203,  the  Lewis  acidity  of  the  A\(Td)  site  is 
better  represented  by  the  Al32048  cluster,  whereas  the  AlgOi2 
overestimates  the  acidity  and,  therefore,  A Edis.  On  the  other 
hand,  the  A180i2  cluster  gives  a  better  description  of  the  Lewis 
acidity  (and  thus  of  AEads)  for  the  Al(7d)  site  on  the  surface  of 
the  nonspinel  model.84  The  unpassivated  A12o03o  cluster23  (i.e., 
with  no  PHT)  significantly  overestimates  the  Al(7)/)  Lewis 
acidity,  and  therefore  the  absolute  A/fads  values.25  However, 
relative  A£ads  values  as  well  as  adsorption  geometries  and  trends 
in  A Z?ads  for  different  bonding  configurations  are  essentially 
independent  of  the  cluster. 

For  DMMP  and  Sarin,  adsorption  via  dative  bond  formation 
between  the  phosphonyl  O  atom  and  the  Al(Td)  is  energetically 
favored  over  other  bonding  modes.  For  example,  adsorption  of 
Sarin  on  the  A180i2  cluster  via  bond  formation  between  Al(7),) 
and  the  F  or  the  isopropoxy  O  atom  gives  A£ads  =  —26.2  and 
—31.0  kcal/mol,  respectively.  These  have  not  been  corrected 
for  BSSE,  and  doing  so  would  decrease  the  magnitude  of  A Z?ads 
by  about  2  kcal/mol.  Both  are  already  much  smaller  than  the 
BSSE-corrected  value  of  —51.5  kcal/mol  (Table  2)  for  adsorp¬ 
tion  via  an  A1  •  •  •  0=P  bond.  The  VX  results  are  discussed  in 
more  detail  below.  Table  2  also  gives  results  for  Av(  P=0 ) 
computed  for  the  smaller  cluster.  A  similar  analysis  for  the  larger 
cluster  is  not  computationally  feasible  at  present.  For  DMMP, 
the  purely  harmonic  values  are  in  good  agreement  with  the 
experimental  result57a,59a  of  —60  cnT1.  For  Sarin  on  the  other 
hand,  the  computed  Av(P=0)  of  —103  cnT1  differs  significantly 
from  the  experimental  value51,62  of  —63  cnT1.  This  discrepancy 
has  been  discussed  previously.25 

The  effect  of  the  adsorbate  on  the  cluster  has  also  been 
considered.  For  the  bare  A180i2  cluster,  with  all  atoms  initially 
in  bulk-lattice  positions,  relaxation  of  the  Al(-0-)3  adsorption 
site  results  in  an  inward  displacement  of  dz  =  —0.26  A  for  the 
A](Td)  which  compares  well  with  the  2-DPS  result76  of  —0.33 
A.  Smaller  displacements  are  found  in  both  cases  for  the  O 
atoms.  Adsorption  largely  reverses  this  inward  displacement. 
After  adsorption  of  DMMP,  for  example,  dz  for  the  A  1(7’,,)  is 
only  —0.04  A  relative  to  the  unrelaxed  bare  cluster.  Thus, 


adsorption  causes  an  outward  displacement  of  about  0.22  A 
relative  to  the  relaxed  bare  cluster. 

For  VX,  bonding  between  A\(Td)  and  the  O  atom  of  the 
CH3CH20-  group  is  not  considered  here  since  this  is  already 
known,  in  the  case  of  DMMP  and  Sarin,  to  be  less  favorable 
than  bonding  to  the  O  atom  of  the  P=0  group.  However,  it  is 
necessary  to  investigate  adsorption  via  an  A 1  •  •  •  S  or  A1  •  •  •  N 
bond.  The  latter,  in  particular,  might  be  strong  at  a  Lewis-acid 
surface  site.  The  use  of  the  small  A180i2  cluster  reduces,  as 
much  as  possible,  the  importance  of  steric  factors  and  of  multiple 
bonding  interactions  (see  below).  This  enables  an  unbiased 
evaluation  of  the  electronic  contribution  of  just  a  single  dative 
bond  to  the  adsorption  energies.  The  results  (not  corrected  for 
BSSE)  are  AEads  =  —30.4  kcal/mol,  r(Al— N)  =  2.10  A  for 
A1***N  and  A£ads  =  —34.8  kcal/mol,  r(Al— S)  =  2.38  A  for 
A1""S.  These  energies  are  significantly  less  than  the  corre¬ 
sponding  BSSE-corrected  value  of  —51.6  kcal/mol  (Table  2) 
value  for  Al  ••  *0=P  formation. 

Figure  5  shows  the  three  optimized  VX  structures.  Passivation 
of  unsaturated  edge  sites  (accomplished  here  through  the  use 
of  PHT)  is  essential  in  the  case  of  a  large,  polyfunctional,  and 
“flexible”  molecule  like  VX.  Results  (not  shown)  for  an 
unpassivated  Al20O30  cluster  indicated  a  tendency  for  the 
molecule  to  distort  so  as  to  allow  an  H  atom  in  one  of  the 
isopropyl  — CH3  groups  to  interact  with  an  Al  edge  site.  With 
PHT,  this  spurious  effect  is  not  seen  even  for  the  smaller  A18Oi2 
cluster.  When  adsorbed  on  either  cluster,  all  three  species 
exhibited  two  or  more  CH  •  •  •  O  bonds  between  alkyl  H  atoms 
and  O  atoms  on  the  A1203  surface  which  were  identified  using 
AIM  analysis.  This  form  of  H-bonding  is  commonly  found  in 
the  adsorption  of  phosphonyl  species  on  oxide  surfaces. 19-22  26 
The  AIM  results  for  the  hydroxylated  surfaces  are  discussed  in 
more  detail  below. 

5.2.  Adsorption  on  Hydroxylated  Surfaces.  5.2.1.  Al80I2 
Cluster.  Experiment57'59,62  suggests  that  Al— OH-"0=P  bonding 
(in  addition  to  A1***0=P  dative  bonding)  may  be  important 
when  OH  groups  are  present.  The  first  set  of  calculations  was 
done  using  the  A180i2  cluster  so  that  Ar(P=0)  could  be 
obtained  and  so  that  bonding  to  different  types  of  OH  groups 
could  be  studied  individually  while  avoiding  as  much  as  possible 
any  complications  due  to  other  forms  of  bonding. 

5. 2. 1.1.  DMMP.  Results  for  the  adsorption  of  DMMP  will 
be  analyzed  in  detail  first,  in  order  to  make  comparisons  with 
experiment  and  to  provide  a  framework  for  subsequent  discus¬ 
sions  of  Sarin  and  VX,  for  which  fewer  experimental  results 
are  available.  The  surface  was  hydroxylated  in  the  1,2  or  1,4 
configuration  (Figure  3a),  and  a  DMMP  molecule  was  then 
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Figure  5.  Optimized  structures  for  bonding  of  VX  to  the  Alg0i2  cluster 
(with  PHT).  Bonding  is  via  (a)  the  N  atom,  (b)  the  S  atom,  and  (c)  the 
phosphonyl  O  atom.  Configuration  (c)  is  energetically  favored  (see  text). 
For  clarity  the  PHs  are  not  shown. 

positioned  with  the  O  atom  of  the  P=0  or  P— O— CH3  group 
oriented  toward  one  OH  or  the  other  and  the  geometry 
optimized.  The  optimization  procedure  was  the  same  as  that 
described  above  for  the  hydroxylated  surface  but  also  included 
the  DMMP.  The  resulting  adsorption  energies  are  given  in  Table 
3,  and  some  representative  structures  are  shown  in  Figure  6. 
Most  of  the  following  discussion  will  focus  on  the  type  Illb 
structure  since,  as  noted  above,  this  is  more  accurately  described 
by  the  cluster  models  than  is  the  type  Ilia. 

For  the  P=0  group,  A£ads  =  —4.5  and  —23.2  kcal/mol  were 
found  for  H-bonding  to  the  type  I  and  type  Illb  site,  respectively, 
with  r(H-"0)  distances  of  1.837  and  1.451  A.  Of  the  two.  the 
type  III  is  clearly  more  favorable  which  is  consistent  with  its 
higher  Brpnsted  acidity  and  with  experimental  results57,59,62 
showing  that  type  I  OH  sites  are  not  the  primary  sites  for  the 
initial  adsorption  of  DMMP  or  Sarin.  However,  even  the  highest 
H-bond  energy  is  smaller  than  A£ads  for  Al(7rf)  •  •  •  0=P  bonding 
to  an  OH-free  cluster  (Table  2).  Thus  the  P=0  group  is  favored 
both  for  H-bonding  and  for  dative-bonding  to  the  Al(7j/). 

Figure  6  shows  that,  in  addition  to  the  main  O— H"‘0=P 
bond,  some  of  the  methyl  H  atoms  form  C— H*"0  bonds.  These 
were  examined  using  the  AIM  approach,  and  some  results  are 
given  in  Table  4.  The  application  of  AIM  specifically  to 
C— H-"0  bonds  has  been  discussed  previously. 90-92  The 
C— H-  •  •  O  bonds  found  here  all  exhibit  3,-1  critical  points39,40 
with  a  charge  density  (p)  and  a  Laplacian  of  the  charge  density 
(V2p)  in  the  range  expected90  for  such  bonds.  The  bond  lengths 
and  angles  are  also  characteristic  of  C— H*"0  bonds  90,93,94  It 
has  been  shown90,95  that  the  bond  energy  between  two  species 
forming  more  than  one  H-bond  is  proportional  to  the  total  p  at 
all  H-bond  critical  points.  On  this  basis,  Table  4  shows  that  the 
O— H"*0  bond  accounts  for  about  70%  of  AEld:,.  The  C—  H •  •  •  O 
bonding  is,  therefore,  not  sufficiently  strong  to  have  a  significant 


effect  on  the  comparison  between  A.Eads  for  different  bonding 
configurations. 

In  order  to  check  for  basis-set  dependence,  some  calculations 
were  repeated  using  TZ2P  basis  sets.  These  differ  from  the  TZP 
basis  sets  described  above  in  the  addition  of  a  second  polariza¬ 
tion  shell  (d  for  H  and  f  for  all  other  elements).  The  calculations 
were  done  for  adsorption  via  an  OH"  •0=P  bond  at  either  the 
type  I  or  type  Illb  site.  The  calculations  include  the  geometry 
optimization  of  the  bare  cluster,  the  free  molecule  and  the 
cluster+molecule  combination  as  well  as  the  BSSE  correction. 
The  results,  —3.7  and  —24.3  kcal/mol  respectively,  are  close 
to  the  corresponding  TZP  values  of  —4.5  and  —23.2  kcal/mol 
given  above. 

Another  check  involved  the  choice  of  functional.  Zhao  and 
Truhlar96  have  designed  a  series  of  GGA  XC  functionals  which, 
for  H-bonding  and  noncovalent  interactions,  are  often  superior 
to  those  more  commonly  used  (such  as  the  PBE  functional 
employed  here).  In  ADF  2007  these  functionals  can  be  applied 
a  posteriori  using  the  SCF  density  obtained  in  a  PBE  calculation. 
The  M06— 2X  functional96  was  used  to  re-evaluate  AEads  for 
DMMP  bonding  via  the  P=0  group  to  either  the  type  I  or  type 
Illb  OH  site  (using  TZP  basis  sets).  The  results  (uncorrected 
for  BSSE)  were  —5.7  and  —29.6  kcal/mol,  respectively,  for 
M06— 2X  vs  —5.1  and  —24.5  kcal/mol  for  PBE  (also  uncor¬ 
rected  for  BSSE).  The  results  are  not  sufficiently  different  to 
affect  the  present  analysis. 

Vibrational  frequencies  were  computed  for  DMMP  adsorbed 
at  the  type  Illb  OH  site.  A  Av(P=0)  of  —78  cm-1  was  found, 
relative  to  the  computed  gas-phase  value.  This  is  a  somewhat 
larger  shift  than  is  seen  for  adsorption  via  an  Al(7);) ••  •  0=P 
bond,  either  in  experiment572,592  (—60  cm-1)  or  computationally 
(—62  cm-1.  Table  2).  Apparently  H-bonding  can  have  a  larger 
effect  than  dative  bonding  on  v(P=0),  even  though  A£ads  is 
smaller.  Bertilsson  et  al.49  have  reported  IR  data  for  DMMP 
adsorbed  from  the  vapor  phase  onto  alkanethiolate  self- 
assembled  monolayers  (SAMs)  terminated  in  —OH  or  — C02H 
tail  groups.  Adsorption  occurs  via  an  OH  •  •  •  0=P  bond  analo¬ 
gous  to  those  studied  here.  For  bonding  to  the  acidic  — C02H 
group,  v(P=0)  was  found  at  1220  vs  1276  cm-1  in  the  gas 
phase,  giving  Au(P=0)  =  —56  cm-1.  A  room-temperature 
adsorption  enthalpy  of  A Ha  15.3  kcal/mol  (64.0  kj/mol)  was 
estimated  by  applying  Brunauer— Emmett— Teller  theory  to 
adsorption  measurements.  Both  Ai'(  P=0)  and  A//a  are  smaller 
in  magnitude  for  bonding  to  a  terminal  (alcoholic)  OH  group. 
A  gas-phase  thermochemical  calculation  performed  here  (see 
below)  gave  A 7/a  =  21.8  kcal/mol  for  adsorption  via  an 
0H"*0=P  bond  to  the  type  Illb  site  (adopting  the  sign 
convention  of  ref  49).  As  noted  above,  ~70%  of  A.Eads  for 
DMMP  adsorption  on  the  hydroxylated  AlgOi2  cluster  is  derived 
from  the  OH  •  •  •  0=P  bond.  Reducing  the  AEads  contribution  to 
A//a  by  30%  gives  A//a  ~  15.3  kcal/mol  which,  perhaps 
fortuitously,  is  identical  to  the  experimental  estimate  for  bonding 
to  the  —  C02H-terminated  SAM. 

The  preceding  discussion  presents  both  theoretical  and 
experimental  evidence  indicating  that  there  is  no  simple,  direct 
correlation  between  AEads  and  Au(P=0)  for  the  Al(7’d)",0=P 
dative  bond.  This  observation  has  been  made  in  a  previous 
study25  where  a  similar  Av(P=0)  was  computed  for  two  species 
(C13P=0  and  (CH3)3P=0)  with  very  different  AZ?ads  values.  This 
appears  to  be  a  general  result,  suggesting  that  Av(P=0)  depends 
on  other  factors  beyond  just  the  bond  energy.  Understanding 
this  effect  will  require  further  investigation. 

The  computed  harmonic  An(P=0)  values  are  not  sufficiently 
accurate  that  H-bonding  can  be  ruled  out  in  the  experimental 
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DMMP  Sarin  VX* 

OH  type  P=0  P-0-Me  P=0  P-0-Me  P -F  P =0  P-0-Et  P-S-R  R'-2V-(R")2 

type  I  -4.5  -5.0  -4.0  -2.4  +0.34c  -3.8 

type  Ilia  -13.7  -9.9 

type  Illb  -23.2  -13.9  -21.9  -10.4  -8.8  -19.0  -8.0  -10.1  (d) 

“  All  energies  are  in  kcal/mol  and  have  been  corrected  for  BSSE.  The  terms  “type  I”,  etc.  are  defined  in  the  text.  Missing  entries  are  judged 
not  to  be  relevant  (see  text)  and  have  not  been  computed.  The  bold,  italicized  atom  is  the  one  to  which  the  H-bond  is  formed.  b  R"  is  a 
— C(H)(CH3)2  group;  R'  is  the  — CH2CH2—  group;  R  is  the  — R'N(R")2  group.  “Et”  refers  to  the  ethyl  group.  See  Figure  lc.  17  Essentially  no 
bonding  occurs  in  this  case.  d  Geometry  optimization  led  to  proton  transfer  from  the  OH  to  the  N  atom  (see  text). 


Figure  6.  Optimized  structures  for  DMMP  adsorbed  on  the  AlgO^ 
cluster,  (a)  Bonding  at  the  type  Illb  site  to  the  P=0  group;  (b)  Bonding 
at  the  type  Illb  site  to  the  P— O— CH3  group;  (c)  like  (a)  but  bonding 
at  the  type  Ilia  site;  (d)  like  (b)  but  bonding  to  the  type  Ilia  site.  The 
numbered  dashed  lines  show  the  H-bonding  interactions,  corresponding 
to  Table  4.  (1)  and  (1')  are  OH---0  bonds  and  (2)  and  (3)  are  CH---0 
bonds.  In  (d),  bond  (2)  is  partially  obscured  by  the  A1  (Tj)  atom. 

data,57a,59a  for  DMMP  adsorption  on  partially  hydroxylated 
y-Al203,  strictly  on  the  basis  of  the  better  agreement  between 
the  observed  Av(P=0)  and  that  computed  for  dative  bonding. 
Nevertheless  the  results  are  consistent  with  the  conclusion,  based 
on  AEads,  that  dative  bonding  is  favored  over  H-bonding  when 
vacant  AlfT)/)  sites  are  available.  Similar  conclusions  were 
reached  by  Templeton  and  Weinberg593  who  found  that  A1*”0=P 
bonding  is  energetically  favored  over  OH  •  •  •  0=P  bonding  but 
gives  about  the  same  Ai'(P=0).  H-bonding,  although  relatively 
weak,  can  occur  when  OH  groups  are  present  and  is  the  only 
mode  of  adsorption  on  a  fully  hydroxylated  surface.  The  IR 
spectrum62  of  Sarin  on  y-Al203  is  essentially  independent  of 
whether  or  not  OH  groups  are  removed  by  evacuation  at  900 
°C  prior  to  adsorption.  This  also  suggests  that  adsorption  can 
occur  at  both  A1  and  OH  sites  and  that  the  two  cannot  easily  be 
distinguished  solely  in  terms  of  the  internal  vibrational  modes 
of  the  adsorbate. 

Finally,  a  calculation  was  performed  for  adsorption  at  the 
type  Illb  site  of  the  Al32048  1,4  model  in  Figure  3b.  H-bonds 
to  the  O  atom  of  both  the  P=0  and  P— O— CH3  groups  were 
investigated.  In  either  case  geometry  optimization  resulted  in 
movement  of  the  molecule  away  from  the  OH  site  and  toward 
one  of  the  vacant  A\(Td)  sites  on  the  periphery  of  the  cluster, 
where  a  dative  bond  then  formed.  This  illustrates  a  point 
discussed  above;  namely,  that  such  sites  are  preferred  for 


adsorption  when  they  are  available  (i.e.,  when  not  all  Al(7j/) 
have  formed  type  I  sites).  The  results  also  indicate  the  reason 
why  this  initial  series  of  calculations,  aimed  mainly  at  a 
comparison  of  adsorption  at  different  types  of  OH  sites,  was 
performed  using  the  small  A180i2  cluster  for  which  the  only 
surface  Al(7j/)  site  is  occupied  by  an  OH. 

5. 2. 1.2.  Sarin  and  VX.  The  adsorption  of  Sarin  and  VX  on 
the  AlgOi2  cluster  was  evaluated  in  the  same  manner  as  for 
DMMP.  In  view  of  the  DMMP  results,  attention  was  focused 
on  interactions  between  the  type  Illb  OH  and  either  type  of  O 
atom  (alkoxy  or  phosphonyl)  and  on  the  bonding  of  F  (Sarin), 
N  (VX),  or  S(VX)  to  either  of  the  two  OH  sites.  The  results, 
shown  in  Table  3,  are  consistent  with  those  for  DMMP.  Bonding 
via  the  interaction  of  the  P=0  group  with  the  type  Illb  OH  is 
the  most  favorable  mechanism.  The  AIM  analyses  (not  shown) 
were  consistent  with  those  in  Table  4.  The  trends  in  p  and  V2p 
at  the  3,-1  bond  critical  points  reflected  those  in  AEads,  with 
the  major  contribution  coming  from  the  bond  between  the  type 
Illb  OH  and  the  O  atom  of  the  phosphonyl  or  the  alkoxy  group. 
In  the  case  of  Sarin  adsorption  via  a  weak  OH  •  •  •  F  bond,  a  large 
fraction  of  the  small  AEads  derived  from  an  additional  CH  •  •  •  O 
bond.  Bonding  of  VX  via  an  OH  •  •  •  S  bond  at  the  type  I  site 
gave  a  very  slightly  endothermic  AEads  of  +0.34  kcal/mol 
suggesting  a  very  weakly  bound  metastable  state.  A  more 
accurate  description,  in  this  case,  may  require  geometry 
optimization  using  a  GGA-XC  functional96  appropriate  to 
noncovalent  interactions. 

It  is  noted  that  the  trend  in  the  type  Illb  OH*--0=P  bond 
energy  (DMMP  >  Sarin  >  VX)  parallels  that  for  the 
Al(7’d)",0=P  bond  on  the  OH-free  surface  (Table  2).  However, 
the  energy  differences  are  in  all  cases  small,  and  small  secondary 
contributions  from  C— H"*0  interactions  make  a  detailed 
comparison  difficult.  Calculations  of  the  vibrational  modes  gave 
Av(P=0)  =  —63  and  —84  cnT1  respectively,  relative  to  the 
computed  gas-phase  frequencies,  for  Sarin  and  VX  adsorbed 
via  a  type  Illb  OH  •  •  •  0=P  bond.  These  are  comparable  to  the 
corresponding  value  of  —78  cm-1  obtained  above  for  DMMP. 
The  VX  result  is  also  close  to  the  value  of  —95  cnV1  (Table  2) 
obtained  for  dative  bonding  to  the  OH-free  surface. 

An  interesting  effect  was  observed  when  attempting  to  adsorb 
VX  via  an  OH  •  •  •  N  bond  at  the  type  Illb  site.  Geometry 
relaxation  led  to  transfer  of  the  proton  to  the  N  to  give  the 
corresponding  ammonium  ion.  The  reaction  energy  of  —18.8 
kcal/mol  (uncorrected  for  BSSE)  was  somewhat  less  than  the 
(uncorrected)  value  of  A£ads  =  —24.5  kcal/mol  for  the  type  Illb 
OH  •  •  •  ()=P  bond.  This  is  an  unexpected  result,  given  that 
experiment97  shows  that  adsorption  of  pyridine  on  hydroxylated 
y-Al203  does  not  yield  a  pyridinium  ion.  This  led  to  an 
examination  of  the  proton  affinities  (PAs)  of  pyridine,  2,6- 
dimethyl  pyridine  (2,6-DMP)  and  VX.  2,6-DMP  is  of  interest 
since  it  gives  evidence97  of  being  a  sufficiently  strong  base  to 
form  the  corresponding  pyridinium  ion  when  adsorbed  on 
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TABLE  4:  AIM  Analysis  of  Selected  H-Bonded  Structures  for  Adsorbed  DMMP 

Bermudez 

structure 

A-^ads 

bond 

r(O---H)  b 

r(C  •  •  •  O)  b 

X(O---H-X)  c 

P 

V2p 

Figure  6 a 

-23.2 

1 

1.451 

166 

0.081 

0.129 

2 

2.143 

3.216 

163 

0.018 

0.061 

3 

2.295 

3.349 

160 

0.013 

0.046 

Figure  6b 

-13.9 

1 

1.580 

170 

0.062 

0.128 

2 

2.351 

3.333 

147 

0.012 

0.041 

3 

2.261 

3.254 

149 

0.015 

0.051 

Figure  6c 

-13.7 

1 

1.493 

174 

0.072 

0.130 

2 

2.161 

3.165 

151 

0.017 

0.062 

3 

2.394 

3.419 

155 

0.011 

0.039 

Figure  6d 

-9.9 

1 

1.670 

174 

0.048 

0.119 

r 

2.105 

148 

0.018 

0.066 

2 

2.340 

3.409 

163 

0.012 

0.042 

3 

2.269 

3.337 

164 

0.015 

0.048 

“  AE^S  in  kcal/mol  (cf.  Table  3);  p  (charge  density  at  the  “3,-1”  critical  point)  and  V2p  (Laplacian  of  charge  density  at  the  critical  point)  in 
atomic  units.  The  bond  numbers  refer  to  Figure  6.  b  Interatomic  distances  in  A.  r(0-”H)  is  the  length  of  the  H-bond.  r(C*"0)  is  the  distance 
from  C  to  O  in  the  C— H---0  bond.  c  0”‘H— O  (bond  1,T)  or  0'"H— C  (bond  2,  3)  angle  in  degrees. 


y-AhO;.  The  PA  of  an  amine  XN  at  finite  temperature  (7)  is 
defined  as  the  enthalpy  change,  AH(T),  for  the  gas-phase 
reaction  XN  +  H+  —  XN— H+  and  is  given  by 

A£pA  =  t(XN-H'  )  -  E(XN)  -  £(H+) 

AH(T)  =  A £PA  +  £trv(XN-H+)  -  £TRV(XN)  -  (5/2 )RT 

Here  F(XN— H+)  and  F(XN)  are  the  relaxed  electronic 
energies  of  the  protonated  and  the  bare  amine,  and  ETRV  is  the 
total  translational,  rotational  and  vibrational  energy.  The  (5/ 
2)RT  accounts  for  the  H+  translational  energy  and  for  the  A(PV) 
contribution.  In  this  definition,  a  negative  AH{T)  indicates  an 
exothermic  process. 

The  experimental98  PAs  of  gas-phase  pyridine  and  2,6-DMP 
at  room  temperature  are  A//(298)  =  —930  and  —963  kJ/mol 
respectively.  The  values  computed  here  for  pyridine,  2,6-DMP 
and  VX  are  —926.0,  —963.2,  and  —986.4  kJ/mol  respectively. 
The  pyridine  result  is  also  in  good  agreement  with  the  value  of 
-935.5  kJ/mol  obtained99  at  the  B3LYP/6-311+G(2df,p)  level. 
The  computed  A/ /( 298 )  for  the  amine  group  in  VX  is  also  close 
to  the  experimental  values98  for  other  trialkyl  amines  such  as 
triethyl-  and  tripropylamine  (—981.8  and  —991.0  kJ/mol 
respectively).  These  results  indicate  that  the  VX  amine  group 
is  more  basic  than  2,6-DMP  and  should  thus  be  able  to 
deprotonate  acidic  sites  on  the  y-AhCh  surface.97  It  is  noted 
that  AZidp  for  the  OH  site  (Table  1)  is  larger  in  magnitude  than 
AEfA.  However,  A Z?dp  as  defined  above  includes  a  large 
contribution  (on  the  order  of  several  hundred  kJ/mol)  from  the 
energy  required  to  separate  the  Al— O  and  the  H+  to  infinity. 
Hence  AF^  greatly  overestimates  the  energetic  cost  to  move 
the  H+  to  the  nearby  (adsorbed)  VX  molecule. 

5.2.2.  Al32048  Cluster.  This  section  builds  on  the  work 
described  above  by  treating  adsorption  on  a  fully  hydroxylated 
surface,  one  in  which  multiple  OH  •  •  •  0=P  and  CH  •  •  •  O  bonds 
are  possible.  Previous  work24  on  the  adsorption  of  phosphonyl 
species  on  the  hydroxylated  a-Si02  surface  indicated  that 
formation  of  two  Si— 0H*"0=P  bonds  to  a  single  P=0  group 
is  the  energetically  favored  mode  of  adsorption.  It  is  of  interest 
to  determine  whether  this  model  can  be  extended  to  y-AI20> 

These  calculations  employed  the  AI32O48  cluster  with  all  seven 
“surface”  M(Td)  sites  (cf.  Figure  3b)  converted  to  type  I  OH 
sites.  This  was  done  to  inhibit  the  preferential  formation  of  an 
Al(7)/)  •  •  •  0=P  bond  which  occurs  during  adsorption  of  DMMP 
when  empty  Al(Frf)  sites  are  available  (see  above).  The  seven 
accompanying  H  atoms  were  placed  on  surface  O  atoms  to  form 
type  III  sites.  Six  of  these  were  placed  on  Oa  sites  to  form  1,2 


configurations  with  the  peripheral  type  I  sites.  The  seventh  was 
placed  on  an  Ob  site  to  form  a  1,4  configuration  with  the  central 
type  I  Al-OH  site.  The  atoms  involved  are  labeled  “(1)”  and 
“Ob(4)”  in  Figure  3b.  This  adsorption  site,  which  is  used  in 
subsequent  calculations,  thus  resembles  that  studied  for  the 
Al80i2  cluster  but  is  now  located  at  the  center  of  a  fully 
hydroxylated  surface.  Geometry  optimization  included  the 
adsorbed  molecule  (if  any),  all  the  added  H  and  OH,  all  seven 
surface  Al  atoms  and  all  surface  O  atoms  to  which  an  H  had 
been  added. 

A  point  of  concern,  for  a  complex  system  such  as  this,  is 
whether  a  true  global  energy  minimum  is  located  in  the 
relaxation  process.  As  a  test,  two  sets  of  DMMP  calculations 
were  done.  One  began  with  a  cluster  which  had  been  previously 
relaxed  when  bare.  The  other  began  with  the  cluster  in  a 
configuration  which  was  far  from  relaxed,  with  all  atoms  in 
bulk-lattice  positions  and  all  O— H  and  Al— O— H  bonds  normal 
to  the  surface.  In  either  case,  bonding  of  the  phosphonyl  O  atom 
to  the  central  type  I  OH  and  to  the  adjacent  type  Illb  site  were 
treated.  For  adsorption  at  the  type  I  site  the  relaxed  energy 
minimum  was  independent  of  the  starting  configuration  to  within 
0.05  kcal/mol.  For  the  type  III  site,  adsorption  on  the  previously 
relaxed  cluster  gave  an  energy  minimum  that  was  slightly  lower 
(i.e.,  more  stable),  by  0.84  out  of  23.2  kcal/mol  (see  below). 

For  DMMP  and  Sarin  the  structures  and  adsorption  energies 
were  very  similar  to  those  obtained  above  for  the  AlgOn  cluster. 
For  DMMP  adsorption  via  an  OH  •  •  •  0=P  bond,  AFads  =  —23.2 
(—7.2)  kcal/mol  was  found  for  the  type  Illb  (type  I)  site.  These 
can  be  compared  to  the  values  of  —23.2  and  —4.5  kcal/mol, 
respectively,  obtained  for  the  hydroxylated  Al80i2  cluster  (Table 
3).  For  Sarin  an  OH  •  •  •  0=P  bond  at  the  type  Illb  site  gave 
AFads  =  —18.2  kcal/mol  vs  —21.9  kcal/mol  (Table  3)  for  the 
Al80i2  cluster.  As  in  the  case  of  the  A18012  cluster,  AIM  analysis 
showed  the  formation  of  a  few  CH  •  •  •  O  bonds  with  ~70%  of 
the  total  AFads  coming  from  the  type  Illb  OH  •  •  •  ()=P  bond. 

The  results  showed  no  indication  of  bonding  of  the  P  =  O 
to  more  than  one  OH  group,  in  contrast  to  the  case  of 
hydroxylated  a-Si02  noted  above.  The  reason  for  this  is  not 
obvious.  It  is  speculated  that,  because  a  type  III  OH  is  back- 
bonded  to  three  Al  atoms,  it  lacks  the  freedom  of  motion 
available  to  an  OH  on  a-Si02,  which  is  back-bonded  to  only 
one  Si.  This  prevents  the  type  III  OH  from  relaxing  into  a 
configuration  in  which  an  optimum  OH  •  •  •  O  distance  and  bond 
angle  can  be  achieved  simultaneously  for  two  H-bonds  to  a 
single  P=0  group. 
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For  the  fully  hydroxylated  AI32O48  cluster,  severe  steric 
hindrance  was  encountered  when  attempting  to  adsorb  VX,  in 
the  vapor-phase  configuration  shown  in  Figure  lc,  via  a  type 
Illb  OH  •  •  •  ()=P  bond.  This  results  from  one  of  the  isopropyl 
groups  preventing  the  phosphonyl  O  atom  from  approaching 
sufficiently  close  to  the  OF!  group  (see  Supporting  Information). 
The  resulting  bonding  interaction  is  negligibly  small  (Alia(js  ~ 
—  1  kcal/mol  before  BSSE  correction),  presumably  because  the 
steric  repulsion  effectively  neutralizes  the  contribution  from  the 
weak  OH- "O  and  C H •  •  •  O  bonds  (identified  in  AIM  analysis) 
that  are  able  to  form.  The  Supporting  Information  discusses 
rotational  isomerization  in  the  isolated  VX  molecule  which,  to 
our  knowledge,  has  not  been  analyzed  previously.  Steric 
hindrance  does  not,  however,  affect  formation  of  the 
Al(7rf)*"0=P  dative  bond  on  the  OH-free  surface  since  in  this 
case  the  bonding  geometry  (cf.  Figure  5c)  keeps  the  amine  group 
well  away  from  the  surface,  even  for  the  larger  AI32O48  cluster. 

In  order  to  facilitate  adsorption  on  the  fully  hydroxylated 
AI32O48  cluster,  the  free  VX  molecule  was  placed  into  a  more 
favorable  configuration,  as  described  in  the  Supporting  Informa¬ 
tion.  This  was  only  ~2  kcal/mol  higher  in  energy  than  that 
shown  in  Figure  lc;  however,  the  configuration  (chosen  to 
minimize  steric  hindrance)  corresponds  to  a  relatively  deep  local 
minimum  in  energy.  A  BSSE-corrected  AE.d(k  of  —11.5  kcal/ 
mol  was  found  for  adsorption  via  a  type  Illb  OH"  •  0=P  bond. 
Unlike  the  case  for  DMMP  and  Sarin  discussed  above,  this  is 
significantly  smaller  in  magnitude  than  the  result  (—19.0  kcal/ 
mol,  Table  3)  for  the  same  adsorption  site  on  the  hydroxylated 
Al80i2  cluster.  Only  a  small  part  of  this  difference  can  be 
explained  by  the  ~2  kcal/mol  difference  in  free-molecule 
energies.  It  is  speculated  that  the  remaining  discrepancy  results 
from  steric  effects,  which  prevent  complete  relaxation  of  the 
adsorbed  molecule,  or  from  the  fact  that  the  adsorbed  VX  may 
be  “trapped”  in  a  local  energy  minimum  lying  significantly 
higher  than  the  true  minimum.  Locating  the  global-minimum- 
energy  for  a  large  molecule  like  VX  on  hydroxylated  y-AF03 
will  require  a  molecular  dynamics  treatment  at  finite  temperature. 

5.3.  Electronic  Excitation.  Tables  SI— S3  in  the  Supporting 
Information  show  TDDFT  results  (transition  energies  and 
oscillator  strengths)  for  dipole-allowed  (singlet  — *  singlet) 
transitions  of  the  gas-phase  species.  The  main  concern  is  with 
transitions  lying  closest  to  the  near-UV  edge  of  the  TSS  (Figure 
7);  hence,  Tables  SI— S3  show  only  the  twelve  lowest-energy 
excitations.  As  noted  above,  there  are  few  experimental  data 
for  the  electronic  excitation  of  these  species.  In  the  absence  of 
such  data,  results  obtained  using  ADF  and  Gaussian  as  described 
above  were  compared.  In  the  latter  case  the  molecular  geometry 
was  first  optimized  at  the  B3LYP/6-31 1++G(d,p)  level.  In 
general  the  agreement  between  the  two  sets  of  TDDFT  results, 
obtained  using  different  functionals  and  basis  sets,  is  quite  good. 
Those  obtained  using  SAOP  are  considered  more  reliable  since 
this  functional  was  specifically  designed44  to  model  electronic 
excitation.  The  results  all  pertain  to  T  =  0  K  with  the  molecule 
frozen  in  the  lowest-energy  ground-state  configuration.  Including 
the  effects  of  finite  temperature100  can  lead  to  shifts  in  the 
transition  energies  which  are  not  taken  into  account  here. 
However,  these  effects  are  not  large  enough  to  preclude  the 
central  goal  of  determining  the  overlap  of  the  computed  spectra 
with  the  TSS. 

For  free  DMMP  the  onset  of  absorption  is  computed  to  lie  at 
7.12  eV.  The  low-energy  transitions  correspond  to  n  — 1 •  a* 
excitations  of  nonbonding  orbitals  (NBOs)  on  the  phosphonyl 
and  methoxy  O  atoms  and  lie  well  outside  the  region  of 
significant  TSS  intensity.  The  results  are  consistent  with  the 
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Figure  7.  The  terrestrial  solar  spectrum  from  the  near-infrared  to  the 
near-ultraviolet  (ref  69)  with  arrows  marking  the  computed  onset  of 
electronic  excitation  (at  0  °K)  for  free  VX  (green),  for  VX  adsorbed 
on  OH-free  3/-AI2O3  (red),  and  for  VX  adsorbed  on  hydroxylated 
y-ALCF  (blue).  The  scale  is  linear  in  energy,  and  the  corresponding 
wavelengths  at  the  end  points  are  given. 

available  data70  which  show  no  absorption  (except  for  possible 
impurities)  below  6.2  eV.  The  Sarin  results  are  similar  except 
that  the  absorption  onset  lies  somewhat  higher,  at  7.38  eV.  For 
free  VX,  on  the  other  hand,  the  computed  absorption  onset  is 
at  4.16  eV.  All  transitions  in  Table  S3  involve  excitation  of 
NBOs  on  the  N  and/or  S  atoms,  and  the  threshold  transition  is 
an  n  — *  ct*  excitation  of  the  N  atom.  For  VX  in  ethanol74  the 
lowest-energy  peak  is  at  4.66  eV.  A  blue  shift  relative  to  the 
gas-phase  is  expected101  in  a  solvent  such  as  ethanol  which  can 
lower  the  energy  of  the  ground-state  by  forming  an  H-bond  with 
the  NBO  on  the  N  atom. 

For  molecules  adsorbed  on  the  OH-free  Al80i2  cluster,  the 
excitation  threshold  was  determined  by  identifying  the  lowest- 
energy  transition  for  which  the  initial  and  final  states  are  MOs 
composed  largely  of  atomic  orbitals  (AOs)  associated  with  the 
molecule  (see  below).  For  DMMP,  the  threshold  shifts  to  7.88 
eV  from  7.12  eV  for  the  free  molecule  (Table  SI).  For  Sarin, 
it  shifts  to  8.13  eV  from  7.38  eV  for  the  free  molecule  (Table 
S2).  The  origin  of  these  shifts  (~0.75  eV  in  either  case)  has 
been  examined  in  terms  of  the  MOs  involved.  For  the  free 
molecule,  the  HOMO  and  HOMO-1  are  nearly  degenerate  and 
are  formed  largely  from  NBOs  on  the  phosphonyl  O  atom.  In 
the  adsorbed  state,  dative-bond  formation  shifts  these  to  lower 
energy  by  1.1  eV  or  more,  and  an  MO  composed  of  NBOs  on 
the  alkoxy  O  atom(s)  then  becomes  the  HOMO.  This  MO,  which 
is  not  directly  involved  in  adsorption,  lies  a  few  tenths  of  an 
eV  lower  in  energy  than  the  free-molecule  HOMO.  The  LUMO 
of  the  adsorbed  molecule,  on  the  other  hand,  lies  slightly  higher 
in  energy  than  that  of  the  free  molecule.  The  net  effect  of  these 
adsorption-induced  shifts  is  the  blue  shift  in  the  excitation 
threshold. 

The  computed  HOMO— LUMO  gap  of  the  bare  Al80i2  cluster 
(with  PHT)  is  4.4  eV  (ref  25)  vs  the  experimental68  y-Al203 
band  gap  of  7.2  eV.  Hence,  the  TDDFT  calculations  yield  many 
low-energy  transitions  between  cluster  levels  as  well  as  charge- 
transfer  transitions  in  which  the  initial  state  is  an  MO  localized 
mainly  on  the  molecule  and  the  final  state  an  MO  localized 
mainly  on  the  cluster  (or  vice  versa).  These  have  been  ignored 
in  analyzing  the  TDDFT  results.  The  MOs  involved  in  the  near- 
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threshold  internal  excitations  of  adsorbed  DMMP  and  Sarin  are 
composed  almost  purely  of  AOs  associated  with  the  molecule 
(as  opposed  to  the  cluster)  with,  in  the  case  of  the  molecular 
LUMO,  a  small  contribution  from  p  orbitals  on  the  AKT-,;) 
adsorption  site.  Thus  the  calculated  threshold  for  internal 
excitation  of  the  adsorbate  should  not  be  strongly  affected  by 
the  DFT  underestimation  of  the  y-Al203  band  gap.  It  is  noted 
in  passing  that  the  excitation  thresholds  for  DMMP  and  Sarin 
are  fairly  close  to  the  actual  y-Al203  band  gap  (7.2  eV). 
Depending  on  the  alignment  of  the  molecular  HOMO  and 
LUMO  with  the  y-Al203  band  edges,  charge-transfer  effects 
may  occur  experimentally  in  photoexcitation  of  the  adsorbed 
species. 

For  VX  (also  adsorbed  via  a  P=0*"  A1  bond)  the  threshold 
red-shifts  to  3.87  eV  from  the  free-VX  value  of  4.16  eV.  The 
threshold  transition  is  the  n  — 1 -  o*  excitation  of  the  NBO  on  the 
N  atom,  as  for  the  free  molecule.  This  orbital,  which  is  not 
directly  involved  in  chemisorption,  is  essentially  unaffected  by 
adsorption.  The  LUMO  of  adsorbed  VX  lies  slightly  lower  than 
that  of  the  free  molecule  which  then  gives  rise  to  the  computed 
red-shift.  For  VX  adsorbed  via  an  OH  •  •  •  Q=P  bond  at  a  type 
Illb  OH  site  the  threshold  is  found  at  4.3 1  eV.  For  completeness, 
excitation  of  a  free  VX  molecule  protonated  on  the  N  atom 
was  also  considered.  As  noted  above,  the  amine  group  in  VX 
is  able  to  undergo  a  Brpnsted  acid— base  reaction  with  an  acidic 
OH  site.  As  expected,  this  eliminates  the  n  — 1 -  o*  excitation  of 
the  N  atom,  and  the  n  —  o*  excitation  of  the  S  atom,  at  5.17 
eV,  then  becomes  the  threshold  transition.  The  VX  results  are 
included  in  Figure  7  which  shows  that  solar  photoexcitation  is 
a  significant  possibility.  It  is  known102  that  n  — *  ct*  excitation 
of  the  N  atom  in  tertiary  amines  leads  to  homolytic  dissociation 
of  alkyl  C— N  bonds  into  free-radical  fragments  which  can  then 
undergo  further  reaction,  and  experiment74  shows  that  irradiation 
of  VX  in  the  vapor  phase  at  4.89  eV  leads  to  molecular 
rearrangement.  The  TDDFT  results  also  indicate  that  protonated 
VX  should  be  significantly  less  photosensitive  to  the  TSS. 

6.  Summary 

The  adsorption  and  photoexcitation  of  DMMP,  Sarin,  and 
VX  on  OH-free  and  hydroxylated  y-Al203  surfaces  have  been 
studied  computationally  using  DFT  and  A180i2  and  Al3204g 
clusters  terminated  with  pseudo-H.  The  design  and  testing  of 
these  clusters  are  important  components  of  the  overall  study. 
The  results  are  as  follows. 

(1)  New  results  for  the  adsorption  of  VX  on  OH-free  y-Al203 
are  consistent  with  previous  results  for  DMMP  and  Sarin. 
Adsorption  via  an  A1  ••  *0=P  dative  bond  at  an  unsaturated 
A1  (TJ)  site  is  favored  over  A1"*N  or  A1-"S  dative  bonding. 
These  are  believed  to  be  the  first  computational  results  for  VX 
adsorption  on  any  surface. 

(2)  Various  issues  have  been  examined  which  are  important 
in  constructing  cluster  models  for  the  hydroxylated  y-Al203 
surface.  These  include  the  configuration  of  different  types  of 
OH  sites  and  the  identification  of  the  sites  involved  in 
chemisorption.  The  models  have  been  evaluated  by  comparing 
ALads  for  dissociative  adsorption  of  H20  with  2-DPS  results, 
and  good  agreement  is  obtained  even  for  the  small  cluster. 
However,  AE^S  shows  a  dependence  on  the  coordination  of  the 
type  III  OH  site  which  is  not  seen  in  the  2-DPS  results. 

(3)  The  valence-band  DOS  for  the  hydroxylated  y-Al203 
surface,  computed  using  a  2-DPS  model,  has  been  compared 
with  photoemission  data.  Type  III  OH  sites  give  a  state  lying 
just  below  the  bottom  of  the  O  2p  band  which  can  be  seen  in 
SXPS.  The  SXPS  feature  was  not  previously  assigned  to  a 


particular  OH  site.  Type  I  sites  give  Al— OH  (AlO— H)  states 
near  the  top  (bottom)  of  the  O  2p  band  which  are  obscured  in 
SXPS  by  the  bulk  A1203  emission.  The  DOS  results  for  both 
the  large  and  small  clusters  are  in  good  agreement  with  those 
of  2-DPS  calculations. 

(4)  Deprotonation  energies  computed  for  the  type  I  and  type 
III  OH  sites  confirm  that  the  latter  is  more  acidic.  Removing  a 
proton  from  the  type  I  site,  followed  by  geometry  relaxation, 
leads  to  reprotonation  via  transfer  of  H+  from  the  type  III  OH 
site.  Similar  deprotonation  energies  are  obtained  for  both  the 
large  and  small  clusters. 

(5)  Bonding  between  a  type  III  OH  and  the  O  atom  of  the 
P=0  group  is  the  favored  mode  of  adsorption  for  all  three 
reagents  on  the  hydroxylated  surface.  Bonding  to  a  type  I  OH 
is  weak,  in  agreement  with  IR  spectroscopy  experiments 
showing  that  type  I  sites  are  not  active  in  the  initial  adsorption 
of  DMMP  or  Sarin.  The  computed  room-temperature  enthalpy 
of  adsorption  of  DMMP  via  a  type  Illb  OH"  •  0=P  bond  agrees 
well  with  an  experimental  estimate  for  DMMP  bonding  to  a 
— C02H- terminated  alkanethiolate  SAM.  However,  even  the 
strongest  OH  •  •  •  ()=P  bond  is  much  weaker  than  an 
AI(7’,/)"'0=P  dative  bond  on  the  OH-free  surface.  The  latter 
mode  of  bonding  will  occur  whenever  vacant  Al(Trf)  sites  are 
available. 

(6)  Unlike  in  the  case  of  hydroxylated  a-Si02  there  is  no 
tendency  of  the  phosphonyl  O  atom  to  bond  to  more  than  one 
OH  group.  It  is  speculated  that  the  type  III  OH  site  lacks 
sufficient  freedom  of  motion  to  permit  the  H-bond  distance  and 
angle  requirements  to  be  satisfied  simultaneously  for  two 
0H*"0=P  bonds  to  a  single  0=P  group. 

(7)  H-bonding  to  the  P=0  group  causes  a  red-shift  in  the 
v(P=0)  stretching  mode.  Ar(P=0)  is  comparable  to  that  seen, 
either  experimentally  or  computationally,  for  Al(7j/),"0=P 
bonding  even  though  the  H-bond  is  considerably  weaker  than 
the  dative  bond.  This  indicates,  in  agreement  with  experiment, 
that  the  two  modes  of  adsorption  cannot  be  distinguished  solely 
by  means  of  the  vibrational  spectrum  of  the  adsorbate. 

(8)  The  amine  group  in  VX  is  able  to  extract  a  proton  from 
an  acidic  type  III  OH  site  to  form  the  corresponding  ammonium 
ion.  The  estimated  reaction  energy  of  —18.8  kcal/mol  is  only 
somewhat  less  than  the  energy  of  the  strongest  H-bond  found 
in  this  study.  Computed  proton  affinities  show  that  the  VX  amine 
group  is  a  stronger  base  than  either  pyridine  or  2,6-dimethyl 
pyridine  and  is  comparable  in  basicity  to  triethyl-  and 
tripropylamine. 

(9)  C— H" •  O  bonds  between  CH3  groups  and  O  atoms,  either 
in  the  A1203  surface  or  at  an  OH  site,  are  formed  in  all  cases 
studied.  However,  AIM  analyses  show  that  these  are  weak  in 
comparison  to  type  III  OH  •  •  •  0=P  bonds. 

(10)  The  energy  surface  of  VX  is  complex,  and  identification 
of  the  global-minimum-energy  configuration  is  difficult  for  either 
the  free  or  adsorbed  species.  Adsorption  on  the  larger  (Al32048) 
hydroxylated  cluster,  and  presumably  on  a  hydroxylated  semi¬ 
infinite  surface,  is  strongly  affected  by  steric  hindrance.  For 
geometric  reasons,  such  factors  are  less  important  for 
Al(7j),"0=P  bonding  on  the  OH-free  surface.  For  DMMP  and 
Sarin  (for  which  steric  effects  appear  to  be  less  important), 
adsorption  on  the  hydroxylated  A180i2  and  Al32048  clusters  give 
similar  results. 

(1  l)For  DMMP  and  Sarin,  either  in  the  gas  phase  or  adsorbed 
on  OH-free  y-Al203,  the  threshold  electronic  excitation  falls  in 
the  vacuum-UV  (~7.1  eV  or  higher),  well  above  the  cutoff  of 
the  terrestrial  solar  spectrum  (~4.5  eV).  However,  VX  excitation 
begins  near  4.2  eV  in  the  gas  phase  and  near  3.9  eV  when 
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adsorbed  via  an  Al(r(/)""0=P  dative  bond  on  the  OH-free 
surface.  Adsorption  via  an  OH  •  •  •  ()=P  bond  shifts  the  VX 
absorption  onset  to  about  4.3  eV,  slightly  higher  than  the  gas- 
phase  value.  Adsorption-induced  changes  in  the  excitation 
thresholds  for  all  three  species  have  been  analyzed  in  terms  of 
shifts  in  MOs.  The  different  excitation  thresholds  for  VX  all 
lie  within  range  of  the  terrestrial  solar  spectrum  which  suggests 
that  photochemistry  induced  by  sunlight  is  a  real  possibility. 
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